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ABSTRACT 

Herein, a detailed protocol for a random mutation 
capture (RMC) assay to measure nuclear point mu- 
tation frequency in mouse tissue is described. This 
protocol is a simplified version of the original 
method developed for human tissue that is easier 
to perform, yet retains a high sensitivity of detec- 
tion. In contrast to assays relying on phenotypic se- 
lection of reporter genes in transgenic mice, the 
RMC assay allows direct detection of mutations in 
endogenous genes in any mouse strain. Measuring 
mutation frequency within an intron of a transcribed 
gene, we show this assay to be highly reproducible. 
We analyzed mutation frequencies from the liver 
tissue of animals with a mutation within the intrinsic 
exonuclease domains of the two major DNA poly- 
merases, S and £. These mice exhibited significantly 
higher mutation frequencies than did wild-type 
animals. A comparison with a previous analysis of 
these genotypes in Big Blue mice revealed the RMC 
assay to be more sensitive than the Big Blue assay 
for this application. As RMC does not require ana- 
lysis of a particular gene, simultaneous analysis 
of mutation frequency at multiple genetic loci 
is feasible. This assay provides a versatile alterna- 
tive to transgenic mouse models for the study of 
mutagenesis in vivo. 

INTRODUCTION 

DNA mutations have been shown to be the underlying 
cause of many diseases including cancer (1,2). Mutations 
can be caused by endogenous oxidative DNA damage, 
damage from genotoxic agents or by error-prone DNA 
replication or repair. Measuring the mutation frequency 



in tissues is essential to understanding the mechanisms 
of mutagenesis in vivo, as well as for measuring the geno- 
toxicity of drugs and compounds. Cellular mutation rates 
can be measured in immortalized fibroblasts (3) or 
lymphocytes (4) by fluctuation analysis. However, this 
method is only amenable to cell types that can be ex- 
planted from the animal, and the mutations being meas- 
ured are generated in vitro in immortalized cells grown on 
plastic, which may behave differently than cells within 
tissues in vivo. 

A sensitive method called random mutation capture 
(RMC) was developed and used to study the frequency 
of random point mutations in tumor and non-tumor 
tissues (5,6). RMC is a genotypic selection technique in 
which genomic DNA is digested with a restriction enzyme, 
followed by single-copy template PCR designed to specif- 
ically amplify molecules with mutations that eliminate a 
particular restriction site. This method allows one to capture 
and selectively sequence rare mutations. The RMC concept is 
based on earlier techniques called Restriction Site Mutation 
(RSM) and restriction fragment length polymorphism- 
polymerase chain reaction (RFLP-PCR) assays (7,8). 

The standard methods in use for measuring the fre- 
quency of random mutations in vivo use transgenic mouse 
models developed for this purpose. The Big Blue (9,10) 
and MutaMouse (11) transgenic models contain lambda 
phage vectors inserted into the mouse genome in a tandem 
array. Mutations in reporter genes within these vectors are 
measured by phenotypic selection after they are shuttled 
out of the mouse genomic DNA, packaged into phage and 
used to infect bacteria. As these assays use engineered 
mouse strains, the transgene must be genetically crossed 
into the mouse genotype(s) of interest before experimen- 
tation. In contrast, the RMC assay is a direct mutation 
detection method that can be done in any mouse model 
without genetic cross or purchase of proprietary mice and 
associated reagents, making it potentially faster and less 
expensive than using transgenic mouse models. 
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Mutations in intron 6 of the human Trp53 locus were 
previously captured using the RMC assay (5,6). Cultured 
human cells had a very low background mutation fre- 
quency, and demonstrated an increase in mutation fre- 
quency with mutagen treatment of the cells that was 
linear with dose (5). The target site demonstrated genetic 
neutrality in that there was no selection for or against 
mutations over time (5). RMC analysis of five human 
tumor and normal tissues indicated a large increase in 
point mutation frequency in tumor as compared to 
normal tissue (6). 

Here, we present a detailed protocol for using RMC to 
measure nuclear mutation frequency in mouse tissue. In 
addition to adapting the method for use in mice, we have 
significantly simplified the original method (5). Removal 
of a lengthy and technically difficult enrichment step 
allows this protocol to be faster than the original 
method, while retaining a high sensitivity. In general, no 
mutations were detectable in young wild-type mice. In one 
mouse where a single mutation was detected, 55 million 
base pairs of DNA were screened giving a mutation fre- 
quency of less than 1 x 10~ 8 mutations/bp. The frequency 
of mutations detected in normal tissues was similar to that 
detected in other endogenous genes, and less than that 
measured by transgenic mouse assays (12). 

DNA polymerase 8 (Pol 8) and DNA polymerase s (Pol 
e) are the two polymerases in the cell that are responsible 
for replication of the bulk of nuclear genomic DNA 
during S phase. These polymerases contain an intrinsic 
3' to 5' exonuclease domain (13). Experiments where this 
domain was mutated have shown that this domain func- 
tions to ensure the fidelity of the replication process in 
yeast (14,15). Homozygous knock-in mice that express 
Pol 8 or Pol s with a point mutation within this exonucle- 
ase domain were created and shown to have higher inci- 
dence of cancer (12,16). In vitro fluctuation analysis 
performed on immortalized mouse embryo fibroblasts 
from these mice demonstrated that mutation of the exo- 
nucelase domains in either of these DNA polymerases gave 
much higher mutations rates (12). In addition, these mu- 
tant mouse stains were crossed into a Big Blue mouse 
background and the mutation frequency within the ell 
transgene was analyzed. In this in vivo mutagenesis assay, 
these mice demonstrated measurable increases in mutation 
frequency in four different tissues (12). Herein, we have 
analyzed these same strains of mice to determine how 
mutation frequency was affected as measured using our 
new RMC assay. Data in the results section illustrate that, 
by RMC, the mice with the polymerase mutant genotypes 
showed a greater increase in mutation frequency over wild 
type than was seen in the Big Blue assay. Additional 
findings based on analysis of the mutation spectra are 
also presented. 

MATERIALS AND METHODS 

Genomic DNA isolation from mouse tissue 

Obtain fresh tissue, or store tissue in vials submerged in 
liquid nitrogen prior to DNA isolation (>100mg best). 
DNA damage is minimized during isolation by not using 



phenol, and by the addition of anti-oxidants and iron che- 
lators to the lysis buffer (17). Use Roche kit for DNA 
isolation (Cat # 11814770001), which uses the high salt 
method for DNA purification. Our procedure follows 
the kit procedure precisely, except for addition of the fol- 
lowing anti-oxidants to the isolation buffer: 4mM histi- 
dine, pH 7.0 (Sigma catalog# 151688; make 100 mM stock 
in water and store at 4°C); 3mM reduced glutathione 
(Sigma catalog# G4705; make 300 mM stock in water 
and store at — 20° C in aliquots) and 1 mM desferrioxamine 
(Sigma catalog# D9533; make fresh 50 mM stock in 
water). Place tissue (fresh or frozen in Nq and quick 
thawed) and lysis buffer in 50 ml screw cap centrifuge 
tube at room temperature. Grind tissue using a tissue 
homogenizer in 5-10 s pulses until all macroscopic pieces 
are gone. Treat with Proteinase K and RNAse as desc- 
ribed in the kit. Transfer to 40 ml polypropylene, high 
speed centrifuge tubes prior to the addition of protein 
precipitation solution and spin to remove insoluble 
protein as described. 

When transferring the DNA containing-supernatant to 
clean 30 ml Corex tube, do not transfer the last bit of 
volume in order to minimize the transfer of any particulate 
material as it will interfere with TaqI digestion later. 
Precipitate DNA with isopropanol and wash pellet with 
70% EtOH as described in kit. Resusupend DNA in 1.4 ml 
TE per 400 mg original wet tissue weight (this will give a 
[DNA] = ~1 mg/ml for liver). Spin samples in a microfuge 
to pellet any remaining insoluble matter and quantify 
using a spectrophotometer. A260/280 should be >1.8. 

Generation of genomic DNA standard curve for DNA 
copy number calculations 

The number of copies per well in the standard curve dilu- 
tions are first approximated by spectra photometric ab- 
sorbance at 260 nm (OD 2 6o)- Measure the OD 2 6o of an 
uncut total genomic DNA sample to obtain a DNA con- 
centration. Genomic DNA rather than plasmid DNA 
should be used for the standard curve. Divide the concen- 
tration of your DNA expressed in pg/ul by the molecular 
weight of the mouse genome (2.74pg/copy) to determine 
the approximate number of copies per microliter in your 
DNA sample. Make your first dilution well above the 
concentration you will use per well in your quantitative 
PCR (qPCR) with experimental primers. We recommend 
starting with 200000 copies/5 ul. Make a 4-fold dilution 
series for 11 dilutions. Mix sample thoroughly before 
making the next dilution. (Note: Make a large volume 
of each dilution — ~ 1 ml per dilution final volume — and 
freeze in single use aliquots.) Store your dilutions at 4°C 
until validated as freeze-thaw cycles alter the effective 
DNA concentration. 

Test your dilution series for linearity in a real-time PCR 
machine with a fluorescence detector (Bio-Rad Opticon 2, 
catalog# CFB-3120 or comparable machine) using the 
control primers (Forward: 5'-CTCACCAAAAACAAAA 
ACAGCC-3', Reverse: 5'-CTTTTGTCCCTCCCACTT 
TGG-3'). For each well in a 96-well plate (Bio-Rad 
catalog# HSP9655), add 12.5 ul SYBR green master 
mix (Applied Biosystems catalog# 4309155), 0.5 ul 
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Uracil-DNA glycosylase (UDG, New England Biolabs 
catalog# M0280L), 0.5 ul of 15 uM stock of control 
primers (Fwd and Rev), 6.5 ul irradiated water and 5 ul 
of diluted DNA. Seal the plate with a plastic seal (Bio-Rad 
catalog# MSB1001), and use the following PCR protocol: 

(i) 37°C lOmin 

(ii) 95°C lOmin 
(hi) 95°C 30 s 

(iv) 60° C 30 s 

(v) 74°C 10 s 

(vi) Read plate. 

(vii) Repeat steps (iii)-(vi) for 47 cycles. 

(viii) Run melting temperature curve to determine relative 
T(m) of products. 

You likely will have to discard a few points from both 
ends of the series, but six or seven dilutions should fall 
on a line in a graph plotting C (T > versus log (copy 
number) as shown in Figure 1. If the series is linear, 
proceed with validation by limiting dilution. 

In order for the absolute copy number to be accurate, 
the approximate DNA copy number values determined by 
optical density are then corrected by limiting dilution 
analysis and estimation of copy number by Poisson distri- 
bution. For limiting dilution analysis, dilutions that ap- 
proximate 1-10 copies per well are tested in multiple wells. 
The number of wells without amplification (those that 
have 0 copies) are used to calculate the average copy 
number per well by Poisson analysis. Select two dilutions 
from your genomic DNA standards (the last one the 
amplified, and the first one that did not amplify — these 
two should be 1-10 copies/5 ul) and run qPCR reactions 
on multiple wells (at least 36 for each of the two dilutions). 
Calculate resultant average copy number per well using 
the Poisson equation: average copy number/well = —In 



[negative wells/total wells]. Multiply average copy 
number/well by a factor of 4 in series to determine absolute 
copy number per well for each DNA standard in the 
curve. See Figure 1 for an example of absolute copy 
number values calculated for a standard curve. Freeze 
your standard curve samples in single use aliquots at 
— 80°C for use on all subsequent experimental samples 
that will be analyzed together. Test one of your frozen 
aliquots before proceeding with control primers and 
mutant-specific primers side-by-side to ensure that both 
primer sets amplify with equal efficiency. The C (T) 
values for a given dilution should be similar for both 
sets of primers with these uncut genomic DNA samples 
(Ideally, within 0.5 C (T) values of each other) for the 
mutation frequency calculation to be accurate. 

TaqI digestion of genomic DNA 

High concentration Taq a I (lOOOOOU/ml; New England 
Biolabs catalog# R0149M) is used to maximize enzyme 
concentration in the reaction. TaqT is a mutant form of 
the enzyme that has a two amino acid replacement at its 
amino terminus. This allows for a higher level of expres- 
sion without interfering with its catalytic properties. 
Achieving sufficient digestion of genomic DNA with 
TaqI is a critical step in this protocol. The quality 
control testing provided by the supplier is designed to 
ensure sufficient digestion efficiency for analysis of 
digested DNA in bulk. For our application, where DNA 
molecules are analyzed individually, we have found that 
the quality of TaqT can vary between lots. Small 
quantities of several lots of TaqT enzyme should be 
tested first with one DNA sample and then the best lot 
then acquired in an adequate amount to complete the ex- 
periment (Figure 2). 
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Figure 1. DNA standard curve for normalization of total genome copy number between tissue samples, (a) A typical display for a qPCR performed 
on standard curve DNA samples run on the Opticon real-time PCR system. In the graph relating C (t) to log quantity, the data points in black are 
the samples that behaved in a linear fashion and were used to generate the standard curve, and those is red were not in the linear range and therefore 
not used, (b) The C (t> values and copy number values (calculated from limiting dilution analysis) are shown in tabular form, color coded to 
correspond to the Opticon graph of amplification curves shown in (a). 
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RMC flowchart 



Test for linearity by qPCR with 
control primers 



Check relative efficiency 

of control and mutant-specific 

primers using DNA standard curve 



IMPORTANT! 

Choose enzyme lot of sufficient 
quality before proceeding 



Order lot of Taq"l that gave the 
most complete digestion 



Ensure near single copy 
template per well so as to 
not mask mutants 



Include only sequence- 
verified mutants in tabulation 



harvest mouse tissues for analysis 



Isolate genomic DNA from tissues 



Generate genomic DNA standard 
curve i 



Test primers and PCR conditions 



Digest a one DNA sample with 
different brands /lots of Taq«l 
restriction enzyme 



Serial dilutions and qPCR to 
test Taq u l cutting efficiency 



Repeat 
If not cut 



Digest the rest of samples with Taq° 
restriction enzyme 



Serial dilutions and qPCR to test 
Taq»l cutting efficiency 




qPCR with mutant- 
specific primers to 
screen for mutants 



qPCR with control 
primers to measure 
total DNA copy number 



Sequence mutant 
products 



Calculate mutation 
frequency 



Figure 2. Flow diagram of RMC protocol. The sequence of experimental steps is shown with important considerations for each of those steps shown 
on the left. Boxed numbers refer to the same numbers step illustrated in Figure 3. A troubleshooting guide for this protocol is available in 
Supplementary Table I. 



TaqI digestion can be done all at once, or over a period 
of several days, storing the samples at 4°C between diges- 
tions. Use one 1.7 ml Eppendorf type tube with 64 ug 
DNA for each sample, bring up the reaction volume to 
800 ul total volume with NEB4 buffer (10X), BSA (100X) 
and double-deionized water. Add 8 ul of Taq a l 
(lOOOOOU/ml) to each tube and digest at 65°C for 1 h, 
shaking at lOOOrpm in an Eppendorf thermomixer or in 
a water bath. Add more Taq"! every 60min for seven 
more hours. Spin samples in a microfuge for 30 s prior 
to each Taq a I addition to remove liquid that collects on 



the lid. A precipitate may form due to denaturation and 
precipitation of the protein in the reaction. With the 
fourth addition or if the samples were held at 4°C over- 
night, give samples a quick microfuge spin, transfer the 
samples to new tubes and add fresh BSA along with the 
Taq a I (to replace denatured, precipitated BSA). (Note: 
loss of BSA carrier protein reduces the efficiency of 
Taq a I digestion.) Next, spin samples in a microfuge and 
transfer to Microcon YM-50 concentrator tubes 
(Millepore catalog# 42416). (Note: do not transfer any 
pellet.) Follow the protocol located in the Microcon 
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mutant-specific primers: 
one plate per sample 
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Figure 3. RMC assay: pictorial illustration of experimental steps. Photograph of a normal mouse liver is shown. Individual steps in the protocol are 
indicated with boxed numbers. In step 2, target sequences in the digested DNA are color coded with blue for wild type and red for mutant. Steps 3, 4 
and 5 illustrate a typical experimental design using 96-well qPCR plates. The dilution factors suggested to test are shown in powers of 4 in the wells. 
In steps 3 and 5, the DNA standard curve is shown in gray, and analysis of six different DNA samples are shown in different colors. In step 4, wells 
marked 'Blk' are blanks to test for DNA contamination during the mutant screening step. Numbered experimental steps shown relate to numbers in 
the experimental flow chart in Figure 2. 



booklet to buffer-exchange the DNA into water, and to 
concentrate the DNA; wash 2 x with 400 ul double- 
deionized water before concentrating the DNA to 
~100 (xl. Bring sample up to 800 ul again in IX Taql 
cutting buffer (i.e. NEB4 plus BSA). Add fresh Taq a I 
enzyme again for two more 60-min incubations at 65°C. 
Perform a buffer-exchange again in the Microcon YM-50 
tubes with 2 x 400 ul washes with UV-irradiated double- 
deionized water, and then concentrate the DNA to a final 
volume of ~100ul. (Note: do not over-concentrate as it 
leads to variable yields.) UV-irradiated, double-deionized 
water is prepared by irradiating 1.0 ml aliquots in open 
Eppendorf tubes in a hood or box equipped with short- 
wave UV light source for 15min. Irradiation destroys any 
contaminating uncut mouse genomic DNA. Finally, dilute 
samples to 800 ul, typically by adding ~700 ul UV- 
irradiated TE (10 mM Tris, pH 8.0; ImM EDTA) and 
store samples at 4°C. (Note: do not freeze DNA samples 
as it alters the effective DNA concentration.) This is con- 
sidered the first 4-fold dilution (4 1 ) as described in the next 
section. 

Testing for Taql digestion efficiency 

Quantitative PCR (qPCR) reactions designed to test the 
efficiency of Taql digestion are performed prior to the 



initiation of screening for mutants. An illustration of a 
typical qPCR plate used to determine Taql cutting effi- 
ciency is shown in Figure 3. Post-PCR restriction analysis 
of products with Taql is used to differentiate between 
wild-type and mutant products. Using samples from the 
Taql digestion protocol above make three 4-fold serial dilu- 
tions with 200 ul DNA and 600 ul IX TE (4 2 , 4 3 , 4 4 ). Mix 
samples thoroughly between dilutions to make the series 
linear. Serial dilutions of the digested DNA are tested with 
mutant-specific primers (Forward: 5'-CTTCTCACCAAAA 
ACAAAAACAGCC-3' , reverse: 5'-ATGGCCAACACACC 
TGACTTCTC-3') by qPCR. [Note: the target Taql site 
used in our protocol in the same one used previously 
(18); however, we used different primers that were 
optimized for this new protocol.] Test triplicate wells 
with each dilution (4 2 ^1 ) for each sample. Run in a 
96-well plate format, using each of the following per 
well: 12.5 ul SYBR Green master mix, 6.75 ul 
UV-irradiated water, 0.5 ul UDG, 0.25 ul of 15 uM stock 
of each primer, along with 5 til of diluted DNA sample. 
Contamination with mouse genomic DNA is always a po- 
tential issue with single copy template PCR so always use 
barrier tips, and change gloves and aliquots of reagent 
frequently, and work, if possible in a hood equipped 
with a short-wave UV light source and irradiate working 
area prior to the start of each experiment. 



e73 Nucleic Acids Research, 2011, Vol. 39, No. 11 



Page 6 of 12 



QPCR products generated with ABI SYBR green mix 
contain some UDP in place of TDP. Addition of UDG to 
the reaction removes any contamination with previously 
generated PCR products. The qPCR protocol to be used 
with the mutant-specific primers is as follows: 

(i) 37°C lOmin (allows for digestion of UTP contain- 
ing DNA). 

(ii) 95°C lOmin 
(hi) 95° C 30 s 

(iv) 62°C 60s 

(v) 74°C 10 s 

(vi) read plate. 

(vii) repeat steps (iii)-(vi) for 47 cycles. 

(viii) run a melting temperature curve to determine 
relative T (m) of products 

When the reaction is done, place qPCR plate immediately 
at — 20°C upon completion as UDG will digest products 
even at 4°C. 

In wells where product is seen, verify that you have the 
correct product by examining the melt temperature curve. 
The correct melting temperature [T (m) ] for products with 
mutant-specific primers is 81.5-82°C. Wells with correct 
T ( m ) product contain either wild-type or mutant product. 
In order to distinguish between wild-type and mutant 
products, perform a post-PCR TaqI digestion and agarose 
gel electrophoresis analysis, using wild-type PCR products 
with and without TaqI digestion as molecular weight 
markers. Digest 8 ul of the each PCR reaction with 
0.25 ul Taq a l (NEB, lOOOOOU/ml), 0.5 ul Uracil-DNA 
glycoslyase inhibitor (UGI, New England Biolabs 
catalog# M0281L), 0.2 ul NEB4 buffer and 1.0 ul 
double-deionized water added per reaction. Incubate for 
lOmin at 65°C, 1 min at 95°C and then 1 min at 25°C 
(digestions are easiest to do in a PCR machine). Analyze 
digested sample using 2% agarose gel electrophoresis. 
Mutant products have one band at 146 bp, while wild- 
type products have two bands at 96 and 50 bp. If the 
products are all wild-type in this test, then your sample 
is not sufficiently digested for analysis. Go back and 
repeat the TaqI digestion step. Otherwise, you will waste 
a lot of reagents and time sequencing wild-type molecules 
that were not cut prior to the mutant screening step. We 
recommend that no more than 10 wild-type molecules 
should be detected per mutant. If mutant products are 
detected, you can use this qPCR analysis of the dilution 
series to determine which dilution will give you approxi- 
mately one mutant copy per well to use for the mutant 
screening (usually 4 2 -4 in our hands). A dilution that 
approximates one copy of template per well is often the 
dilution where no amplification had occurred; amplifica- 
tion will only occur 37% of the time based on a Poisson 
distribution. As a rule of thumb, choose the last dilution 
where some amplification occurred. 

Running qPCR plates to screen for mutants 

Set-up and analyze one full plate for each DNA sample 
using the DNA dilution determined in qPCR plate for 
TaqI cutting efficiency test described above. Use the 



same PCR protocol with the mutant-specific primers 
described above. 

Include 94 wells, all at the same DNA dilution and two 
blank wells with water in place of the DNA template to 
test for DNA contamination. If blanks are contaminated, 
data obtained with this plate must be discarded and 
sample re-run with new reagents. Freeze plates immediate- 
ly at — 20°C after the run as UDG can function quite well 
at 4°C and will degrade your products. 

Analyze the plate for potential mutants. In all wells 
where product is seen, verify that you have the correct 
product by examining the melt temperature curve. 
QPCR analysis with too high of a template copy number 
can cause the masking of mutations with either wild-type 
product or other spurious products if the PCR is not opti- 
mized. To avoid masking mutations, no more that 59 wells 
in each 96-well plate should be positive for product to 
insure near single copy template in every well. 
Determine the sequence of material from all wells with a 
product of the correct melting temperature to identify mu- 
tations. In our experience, given the low rate at which 
mutant molecules are present in the excess of wild-type 
sequences in genomic DNA, a significant number of the 
PCR products obtained will be wild-type in spite of exten- 
sive TaqI digestion. 

Unpurified PCR products can be directly sequenced in 
most high-throughput sequencing facilities in which case, 
2ul or each positive well can be directly sequenced with 
either the forward or reverse PCR primer. Add 0.5 ul 
of UGI to each sequencing reaction to inhibit residual 
UDG activity in sample which will degrade the template. 
In the case where low-throughput or costly sequencing 
is the only option, screening samples by post-PCR TaqI 
digestion and gel analysis as described above can reduce 
the number of samples to be sequenced. 

Determining total DNA copy number 

Set up a 'total copy number plate' using the control primer 
pair and the same PCR protocol described above for the 
generation of the DNA standard curve. Run each plate 
with a standard curve to allow you to convert the C (T) 
values into DNA copy number per well (Figure 1). Up to 
14 DNA samples can be analyzed on one plate for deter- 
mining the total copy number. An example of a total copy 
number calculation plate is shown in Figure 3, step 5. Run 
six wells of each DNA sample, using the same dilution 
that was used in the mutant screening plate. Take the 
average DNA copy number calculated for the six-well rep- 
licates and use this as your total DNA copy number/well 
for that sample. Calculate mutation frequency: mutation 
frequency = (number of sequence verified mutants) / 
(total copies per well) x (number of wells screened) x 
(4 bp per copy). 

Given the stochastic nature of mutagenesis as well as 
the very low frequency of mutation, multiple mutations 
need to be detected before a difference between cohorts 
can be considered real. If two sets of DNA samples being 
compared were to have, for example, similarly low freq- 
uencies as our young wild-type mice presented in the re- 
sults section, more base pairs would need to be screened to 
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determine if there was a difference between cohorts. As a 
rule of thumb, if only 0 or 1 mutation were detected per 
mouse in both sets of samples to be compared, more base 
pairs must be screened until multiple mutants are obtained 
for at least one cohort. A Mann-Whitney statistical test 
can be used to assess statistical significance. 

RESULTS 

Description of the method 

A TaqI restriction site located in intron IV of the Trp53 
locus on chromosome 1 1 was used to develop the mouse 
RMC assay. The target site is illustrated in Figure 4. This 
site is located in an intronic region of no known function 
that is poorly conserved across species, suggesting that it is 
likely to be genetically neutral. Two pairs of primers from 
the Trp53 Intron IV target sequence were selected for use 
that have equivalent amplification efficiencies, one pair 
adjacent to the target TaqI site (control primers) and 
one pair flanking the TaqI site (mutant-specific primers). 
Total genomic DNA is extensively digested with TaqI re- 
striction enzyme. The frequency of uncut wild-type TaqI 
sites that remain in the genomic DNA after digestion must 
be very low in order for rare mutations to be detectable. If 
the frequency of uncleaved, wild-type products captured is 
not in the same range as the frequency of mutant products, 
the vast majority of molecules captured and sequenced 
will be wild type making detection of mutants laborious. 

TaqI digested genomic DNA is screened by qPCR with 
mutant-specific primers to amplify mutant molecules. In 
order to be able to sequence mutant products, the mutant 
screening must be done at high dilution of the DNA 
template, such that only one or two template copies are 
present per well. An illustration of a mutant-screening 
plate is shown in Figure 3. Taql-resistant products are 
verified as mutants by direct sequencing of PCR 
products. The total number of genomes that was screened 
for mutants is measured in parallel for each DNA sample. 
QPCR analysis is performed using the control primers, 
which do not flank the TaqI site, on the same DNA 
dilution that was used for the mutant screening. An illus- 
tration of a typical total copy number plate template is 
shown in Figure 3. As control primers will amplify all 

Mouse Trp53 gene 

KM 011640*1 ■ 1 | | I I | || f ■ MP 03577041 



Fwd M 

FwdC Ta< * 



Pvu II RevC Rev M Rsa I 

Figure 4. Target sequence in Trp53 locus used for RMC assay. A sche- 
matic drawing of the mouse Trp53 locus with the location of RMC 
target sequence is shown. Location of the TaqI restriction site is 
indicated, and control forward (Fwd C), control reverse (Rev C), 
mutant-specific forward (Fwd M) and mutant-specific reverse (Rev 
M) primers are designated with small arrows. 



genomic copies present in each well, the number of total 
DNA molecules screened per well is determined. A DNA 
standard curve is used to convert cycle threshold [C (t>] 
values to DNA copy number. An example of a standard 
curve used for this conversion is shown in Figure 1 . As the 
accuracy of the final mutation frequency number is de- 
pendent on the accuracy of the total DNA copy 
number, it is recommended that a significant effort be 
made to construct a robust set of DNA standards to be 
used on all plates throughout the experiment. 

The final result is a mutation frequency expressed as 
mutations/bp. This number is calculated as: (number of 
sequence verified mutants)/(total copy number per 
well) x (number of wells screened) x (4 bp) 

Reproducibility of RMC analyses 

We analyzed a cohort of 9-month-old, wild-type C57BL/6 
mice using the RMC assay. Mice were sacrificed, and the 
left lobe of each mouse liver was harvested and quickly 
frozen. Genomic DNA was isolated from tissue pieces 
between 200 and 400 mg in weight. In independent experi- 
ments performed 6 months apart, the same DNA samples 
were digested with two different lots of Taq a I enzyme, and 
analyzed with different batches of qPCR reagents. 
Roughly one-fourth of a liver was homogenized to gener- 
ate each DNA sample, and about 1/500 of that DNA 
sample was analyzed in each of the two independent ex- 
periments. Figure 5 shows a comparison of the mutation 
frequencies calculated for the mice from the two different 
experiments. Each bar represents the mutation frequency 
calculated for the individual mice. 

It has been suggested that, as the mutant screening in 
the RMC assay is performed at limiting dilution, signifi- 
cant sampling error is introduced (19). We observed that 
the frequency measured for a given mouse did vary 
between the experiments indicating some sampling error, 
yet the relative difference between the mice was reprodu- 
cible (Figure 5). 

The median mutation frequency for these mice was 
5.3 x 10~ 7 mutations/bp, with >10-fold variation in fre- 
quency between individual mice. These 9-month-old mice 




Figure 5. The reproducibility of data from two RMC analyses of 
9-month-old wild-type C57B1/6 mice. Results from two different experi- 
ments performed on the same DNA samples are presented. Results 
from experiment 1 are shown in black, and those from experiment 2 
are shown in gray. Each bar represents mutation frequencies calculated 
for individual mice. Mouse 'code' names are indicated below each pair 
of data bars. 
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Figure 6. Relative mutation frequency for individual mice from Po/e e,e , Poldl e ^ e and wild-type C57BL/6 genotypes. Mice with different genotypes are 
shown in different colors; Pole e ^ e are black, Poldf^ e in gray and wild type in light gray. Each bar represents data for an individual mouse, which was 
assigned a 'code' at necropsy to blind the experimenter to genotype during analysis. Mouse names are indicated below each datum bar, along with 
the mutation frequency, the number of sequence-verified mutations identified and the total number of base pairs (in millions) screened for that 
individual mouse. 



had significantly higher mutation frequencies than did 
other C57BL/6 mice at 2 months of age shown in 
Figure 6 and described in the following section. These 
results suggest that mice accumulate random mutations 
with age as has been shown previously using other 
methods (20,21). Based on these results, we recommend 
that at least six mice be tested per cohort and that the 
cohorts be age matched. 

Point mutation frequencies in Pol <5 and Pol £ 
exonuclease-deficient mice 

To directly compare the mouse RMC assay to an estab- 
lished method, we analyzed DNA from mice with proof- 
reading deficient DNA polymerase s (Pole e/e ) and DNA 
polymerase 8 (Poldl c/c ). Previously, these mice were 
crossed with C57BL/6 Big Blue mice and analyzed in a 
Big Blue assay (12). Elevated mutation frequencies were 
measured within the ell transgene for both Pole e/e and 
Polclf 1 " Big Blue mice as compared with wild-type 
littermates (12). 

Cohorts of six to seven mice, 6-8 weeks of age, were 
sacrificed from Pole" 1 " and Poldl e/e colonies along with 
wild-type littermates. DNA was isolated from the left 
lobe of the liver, digested with Taq a I and diluted for 
qPCR analysis. For each sample, greater than or equal 
to four million base pairs were screened with 



mutant-specific primers. Results of this analysis are 
shown in Figure 6. Each bar represents the mutation fre- 
quency calculated for the individual mice. Multiple muta- 
tions were detected in all the Pole e/e mice and in six out of 
seven Poldl L ' /e mice. These high mutation frequencies in 
the mutant mice were in sharp contrast to wild-type mice 
where one mutation was detected in only a single individ- 
ual. Thus, elimination of proofreading activity from either 
of the major DNA polymerases caused a significant 
increase in mutation frequency within the RMC target 
TaqI site (Figures 6 and 7). The majority of the 
wild-type mice had no mutations identified giving them 
mutation frequencies of zero. We included less-than 
values shown at the bottom of Figure 6 which reflect the 
number of total base pairs screened for each mouse. In 
order to obtain a background mutation frequency value 
for a wild-type mouse, 55 million base pairs were screened 
from one sample until a mutation was identified. The 
calculated mutation frequency for that one wild-type 
mouse was 9 x 10~ 9 mutations/bp. 

The median mutation frequency in the Pole e/e mice was 
significantly higher than in Poldl" 1 " mice, and both were 
significantly higher than wild type (Figure 7a). Mutation 
frequencies do not have a normal distribution and there- 
fore a standard /-test for parametric data would not be 
applicable. The statistical test we used to analyze these 
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Figure 7. Comparison of RMC analysis of mutation frequency in Trp53 intron IV with analysis of frequency in the ell transgene using Big Blue 
mice, (a) RMC data from the 6- to 8-week-old mice with the different genotypes are shown. Frequencies for individual mice are presented as points 
around the median value that is indicated as a horizontal line. Median values are as follows: 0.0 for wild type, 29.8 for Pole" 1 " and 2.1 for Poldl e/e . 
The median mutation frequencies were found to be statistically different between both mutator mouse genotypes and wild type: wild type versus 
Pole e,e , P = 0.003; wild type versus Poldl ej \ P = 0.009 and Pole' 1 ' versus Poldl el \ P = 0.001. (b) Big Blue mouse data from a previous study (12) on 
6- to 8-week-old mice of the different genotypes. Mutation frequencies measured in four different tissues are shown: small intestine (squares), lung 
(diamonds), thymus (triangles), bone marrow (upside down triangles). The three bars on the left show the effect of a mutation in DNA polymerase 
epsilon, and the three bars on the left show the effects of a mutation in DNA polymerase delta. 



data was a Mann-Whitney test with a correction for ties. 
By performing the 'asymptotic' test rather than the 'exact' 
test, the calculation accounts for tied values. We chose to 
use this test given the large number of zero values 
obtained for the wild-type cohort of mice. Analysis was 
performed using the SPSS 18 statistical program (IBM). 
In Figure 7b, we present a reprinted version of Figure 2A 
from Albertson et al. (12), where they measured the 
mutation frequency in the phage ell gene from Pole e/e 
and Poldl e/e mice using the Big Blue mouse assay. In this 
side-by-side comparison, we see that the relative muta- 
tion frequencies measured in the liver by RMC analysis 
(Figure 7a) were very similar to those measured in four dif- 
ferent tissues using the Big Blue mouse assay (Figure 7b). 
Pole e/e mice had the highest frequency, Poldl c ' le mice had a 
lower frequency, and wild-type mice had the lowest. 

When the median values for the Pole e/e and Poldl 1 '^' 
mice are divided by 249 (the number of mutable base 
pairs in the ell gene) to express them as mutations/bp, 
the values measured by RMC and by the Big Blue assay 
are remarkably similar (27 x 10~ 7 and 26 x 10~ 7 for 
Pole e/e , and 2.5 x 10" 7 and 6.8 x 10" 7 for Poldl e/e , respect- 
ively). However, when the same analysis is done to 
compare the values for the wild-type mice as measured 
by the two methods, RMC gave a lower mutation fre- 
quency than did the Big Blue assay (0.09 x 10~ 7 for the 
one mouse where 55 million base pairs were screened 
versus 2.6 x 10~ 7 ). The larger difference between 
wild-type and mutator mice that we measured with 
RMC as compared to the Big Blue assay may reflect the 
fact that spontaneous mutation frequencies are lower 
within endogenous, transcribed genes. In another study 
where mutation frequencies at the HPRT1 locus and the 
ell reporter gene were compared, a lower frequency of 
spontaneous mutation was measured at HPRT1 than at 



ell in untreated cells, while the mutagen-treated cells gave 
similar values at the two loci (22). 

Comparison of mutation spectra from Pol d and Pol £ 
exonuclease-deficient mice and older wild-type mice 

The very low mutation frequency in young wild-type mice 
precluded our ability to compare the spectra of spontan- 
eous mutations in the polymerase mutant mice with 
wild-type mice from this same study. We therefore used 
the spectra from our RMC analysis of the older cohort of 
wild-type C57BL/6 mice that were used in our reproduci- 
bility study to compare with the polymerase mutant mice 
spectra. In the wild-type mice, a variety of mutational 
events were seen, with the predominant events being 
C-to-T and G-to-A transitions (Figure 8). These are the 
most common mutational events seen in the spectra of 
spontaneous mutations within the ell transgene from the 
Big Blue mouse model (23). It is also the most common 
event seen by deep sequence analysis of human tumor 
tissue (24). The fact that we obtained similar spectra to 
those obtained by other methods strongly suggest that we 
are detecting mutations that occurred in vivo. We detected 
very few T-to-A transversions, which are a common 
mistake made by TaqI polymerase during PCR, suggest- 
ing that a PCR artifact did not contribute significantly to 
our detected mutations. 

C-to-T and G-to-A transitions, which could be con- 
sidered as the same type of mutations occurring on 
opposite DNA strands within the palindromic TaqI site, 
were measured with similar frequency in wild-type mice. 
Poldr /c ' mice had similar spectral profiles as that from the 
wild-type cohort, with the predominant mutational events 
being C-to-T and G-to-A transitions (Figure 8). In 
contrast, the Pole e ' e strain had very large number of 
C-to-T events, but only a single G-to-A event in one 
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Figure 8. Comparison of the spectra of total mutations identified from 
mice with different genotypes. Along the .Y-axis are shown the different 
mutations detected within the TaqI site; each set of stacked bars cor- 
responds to one particular type of mutation detected. The events that 
occurred within a particular mouse are color coded, with one color per 
mouse. 'Code' names for the mice in the legends correspond to those in 
Figures 5 and 6. The spectra shown for the wild-type mice are from a 
9-month-old cohort of C57B1/6 mice used to show the reproducibility 
of RMC analysis shown in Figure 5. These data are included in this 
figure for comparison with the mutator mouse spectra. 



mouse. These results suggest that Pole'' 1 '' mice may have a 
stand bias in the mutational events generated. A similar 
phenomenon was observed with the analogous exonucle- 
ase mutant in yeast (25). Another interpretation of these 
data is that the 'C residue within our target TaqI site is a 
mutational hot-spot for this particular polymerase mutant- 
strain. These data illustrates the potential for RMC, even 
using only one genomic site, to shed light on mutational 
or DNA replication and repair mechanisms. 

Sequence analysis of mutations: independent events versus 
total events 

There are two ways that mutation spectra can be 
tabulated: (i) as 'total' mutations where all mutations 
are included in the tabulation or (ii) as 'independent' 



mutations where, if the same mutation is obtained 
multiple times from the same mouse, it is only counted 
once. The spectra of the 'total' mutations detected by 
RMC from wild-type and polymerase 
proofreading-deficient mice are shown in Figure 8, and 
the spectra for 'independent' mutation can be seen in 
Supplementary Figure 1. For mutagenesis assays with 
larger genetic targets such as those used in transgenic 
mouse assays, duplicate mutations from the same mouse 
are often excluded (26). Given the small size of the TaqI 
target sequence, if only known independent events (differ- 
ent mutations) were counted, it is likely that the number of 
independent mutations would be underestimated. Using 
total or independent events to calculate mutation fre- 
quency both have possible issues. It is best to consider 
the particular experiment and determine the best choice. 
For example, with measuring mutations after acute 
mutagen exposure, one would always count total muta- 
tions as little or no cell proliferation occurred during the 
experiment. With experiments testing for possible mutator 
genotypes that could cause changes in mutation frequency 
throughout the lifetime of the animal, using total muta- 
tions to score for frequency runs the risk of including 
some clonal expansions. In our analysis of the data pre- 
sented here, we chose to count the total mutations rather 
than risk excluding mutations that are repeated simply 
because they are a common event. In our data, the only 
events that were identified more that once from the same 
mouse were C-to-T and G-to-A transitions, which are the 
most common event overall. 



DISCUSSION 

Spontaneous and induced mutagenesis in mice is moni- 
tored in vivo using transgenic mouse models (26,27) or 
in vitro with fluctuation assays using immortalized 
mouse cell lines (9,11,28). The mouse RMC assay 
provides a useful alternative to these other assays for 
measuring mutation frequency, both in vivo and in vitro. 
RMC can be used to study changes in mutation frequency 
as a function of genotype, age or disease state. In addition, 
one could compare the mutation frequencies between dif- 
ferent tissues. 

In transgenic mouse assays, mutations are measured 
within exogenous genes that are transcriptionally silent 
in the mouse and thus would not be subject to 
transcription-coupled repair. Transcription coupled 
repair is a major pathway for cellular DNA repair, espe- 
cially in quiescent cells (29). Comparisons of mutational 
events at the endogenous Hprtl locus and the lad and ell 
transgenes in Big Blue mice revealed that both transgenes 
had a higher spontaneous mutation frequency than did 
Hprtl (22). As RMC does not rely on a reporter gene 
for mutation detection, mutation frequency can be mea- 
sured in any locus in the genome. Using the target TaqI 
site in TRP53 intron IV described in this protocol, the 
mutation frequency within an endogenous, transcribed 
gene is measured and should reflect the rate of mutational 
events in actively transcribed chromatin. Another advan- 
tage of RMC over other random mutagenesis assays is 
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that there is no phenotypic selection step required for the 
identification of mutations. The spectrum of mutations 
sequenced will accurately reflect the events that occurred 
in vivo, not just those that generate a phenotypic change. 

While the Taql target sequence (TCGA) is useful for 
detecting all types of point mutations in that one of each 
the four bases is represented, the size of the genetic target 
in RMC is significantly smaller than that used in the trans- 
genic mouse assays. A Taql site has 4 bp as compared to 
249 bp for the smallest target (ell) used in transgenic 
mouse assays. A smaller genetic target would be predicted 
to decrease the sensitivity of mutation detection. This can 
have the positive effect of lowering the background fre- 
quency in a case where effects of a mutagen or mutator 
genotype are being tested, but could also reduce the ability 
to detect changes at very low mutation frequency. In 
Figure 7, where data from RMC are directly compared 
with those from the ell reporter gene measured in a Big 
Blue mouse assay, a greater difference was measured 
between wild-type and mutator mice with RMC than in 
the Big Blue mouse assay. Thus for this application, RMC 
appeared to be the more sensitive method for detecting a 
mutator phenotype. 

One disadvantage of this RMC protocol as we have 
designed it is that large deletions or insertions are not 
detected. The RMC assay is designed to specifically 
detect point mutations. Given the small size of the PCR 
amplified target DNA, only very small indels are detect- 
able [Figure 8 and see reference (5)]. However, the assay 
can be modified to detect larger indels (30). 

While we describe our protocol for measuring mutation 
frequency in vivo, mouse RMC could also be used to 
measure mutation rate in vitro using immortalized cell 
lines and fluctuation analysis. RMC used in vitro could 
obviate the need for the time consuming drug selection 
step where individual clones of cells must be expanded 
into individual colonies. Fluctuation analysis could be 
performed directly on polyclonal cultures of cells that 
are screened in parallel for random mutations at any en- 
dogenous genetic locus. 

Given that the RMC assay does not require analysis of 
a particular reporter gene, multiple genetic loci could be 
examined in parallel to study the position effect on 
random mutation frequency at the different genomic 
sites. RMC could be used, for example, to compare the 
mutation frequencies between introns and exons, or 
between loci that are transcriptionally active or silent. 
One could also compare euchromatic DNA with hetero- 
chromatic DNA. We describe this protocol using a par- 
ticular genomic site, yet this assay can be adapted to 
examine any Taql restriction site (TCGA), or even 
multiple loci in parallel from the same digested DNA 
samples using different PCR primers. For expanding 
RMC analysis to new genomic sites, a description of im- 
portant design parameters is presented in SI supplemen- 
tary methods. 

Currently, in vivo genotoxic screening of novel drugs or 
chemicals in mice is done using indirect measurements 
such as unscheduled DNA synthesis in the liver (27). 
Tests using either selection or screening for mutations in 
endogenous genes such as HPRT or Dbl are also useful yet 



they are limited to a few tissues. Transgenic mouse assays 
have the advantage of mutation detection in virtually any 
tissue (27), yet their utility is limited due to a high spon- 
taneous mutation rate as compared with endogenous 
genes (22,26). RMC could provide a faster and more sen- 
sitive screening method for both acute and chronic muta- 
genic effects of novel compounds in mice. This simplified 
RMC protocol is amenable to automation, operational in 
any mouse strain, and could also be adapted for use in 
other animal model systems such as rat, dog, or primate. 
Currently, DNA isolation techniques as well as qPCR 
amplification and sequencing have both been success- 
fully adapted to automation. A procedure that combines 
both, coupled with repeated incubations with restriction 
enzyme could allow for an automated RMC assay, which 
would be, by far, the most high-throughput method for 
measuring mutation frequency from cells or tissues. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR online. 
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